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mass is calculated without using the lens equation and compared to a recent observational 
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1 Introduction 



Time delay is one of the four classical tests of general relativity. The first experimental 
confirmation is due to I. I. Shapiro. With his team, he measured in 1968 a time delay of 
240 /is for a (one-way) travel-time of 10 minutes between Mercury and Earth. In October 
2007 Fohlmeister et al. [Ij published a lower bound of 5.7 years on a time delay of truly 
cosmological nature: the travel-time between the quasar SDSS J1004-I-4112 and Earth is 
roughly 10^° years, the lens is a cluster of about 5 ■ 10^'^ solar masses. One month earlier 
Rindler & Ishak [2] corrected the wide held error that the cosmological constant does not 
change the deflection angles of light. Sereno [3] backed up this claim and also derived a 
formula for the time delay with cosmological constant. Earlier results on this time delay 
are due to [U E]. Note also the analysis by Bakola et al. [B] on extreme lensing by black 
holes including a positive cosmological constant. The aim of this paper is to combine the 
two news and to compute the time delay between the best aligned images of the quasar 
(C and D) as a function of the cosmological constant. We assume the idealisation that the 
cluster is static and spherically symmetric and that it dominates all other lensing masses 
in the universe. The fact that there are five images of the quasar of course tells us that 
the cluster is not spherical. 



2 The set up 

Consider a universe, which is empty except for one static, spherical, non-rotating mass 
M, the lens L. A source S*, at rest with respect to the lens, emits photons, which are 
observed at nostra terra T also assumed at rest. We neglect the masses of the source 
and of our local group. We use polar coordinates (r, 6, ip) centered at the lens. Because 
of spherical symmetry the photons' trajectory is in a plane that we take to be ^ = 7r/2. 
The angle ip is measured with respect to the axis defined by the source, ips = 0. With a 
cosmological constant, the gravitational field outside the mass M is given by the Kottler 
metric, 

dr^ = Bdt^- B-'dr^ - rM(^^ 6 = 7r/2, B = l- - lAr^. (1) 

The source has polar coordinates {rs, 0), the earth is at (r^, (Pt). Let a and a' be the two 
physically measured angles between the images and the cluster center and denote by tq 
and Tq the peri-lenses of the two light rays. We write ts and tr for the coordinate times 
of flight from the source to the peri-lens and from the peri-lens to Earth, see figure 1. 
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Figure 1: A double image 
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We will suppose that < 9/10 to avoid the coordinate singularity at the equator 

of the de Sitter sphere. We will also suppose that 6 := GM/tq -C 1, a -C 1 and likewise 
for their primed quantities and keep only terms linear in these four quantities. Let us 
anticipate that in this approximation, a ~ a/1 — Ar'^/Sro/rx- We will also assume that 
A := A/A/3ro ^ 1 and only keep terms linear in this and its primed quantity. For the 
example of the quasar SDSS J1004+4112, all six quantities are of the order of 10~^. 

Our aim is to compute the proper time delay At = ^/bJtt) [t'j, + t'g — — ts) as a 
function of M, a, a', r^, r^, and A. 



3 Integrating the geodesies 



We start with the list of the non-vanishing Christoffel symbols for the Kottler metric with 
9 = n/2 and denote ' := d/dr, 



tr 

r 



B'/{2B), V'u = BB'/2, V\r = - B' / {2B) , 



~rB, r^,,^ = 1/r. 



The geodesic equations read: 



i+B'/Bir = 0, 

f + ^BB'P - iBB'r^ - rBip^ = 0, 
(f + 2r^Vv9 = 0, 



(2) 
(3) 



(4) 
(5) 
(6) 



where we denote the affine parameter by p and ' := d/dp. We immediately get three first 
integrals: 



t 





1/B{r), 
r'^^/Biro) 



Tq B{r) 



B{rQ) 



1/2 



(7) 



(9) 



Eliminating the affine parameter we get: 



dip 
dr 

dt_ 
dr 



± 



1 - 



2GM 2GM r 



r^r'^jrl - 1 L 



-1/2 



To r + To 

2GM 2GM r 



Tq r + ro 



B{r)y^l-r^/r^ 

Integrating equation (ITOi) we obtain a ~ -\/l — Ar'^/Sro/rx and [7] 



-1/2 



(10) 

(11) 



rs aa'Tx 



;i-A4/3)-i. 



(12) 
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Let us integrate equation (|TT|) : 



1 



B{r)^l-rl/r'^ 



1 



2^71/ 2GM r 



ro r + ro 



-1/2 



dr 



~ ro A/-B(ro) [In + 5 /t2 + <5 -^rs + 2 5 /r4] 
We have set ?/ := ro/r, := To/rT > A, 



Its 



dy 



1 



y dy 



2 AVr^A2 

1 



: arctanh 



A / 



Vr^A2 er 



1 



y2 _ A2 /iZT^ 



arctanh 1 



1-A2^ 



(13) 

(14) 
(15) 



1 



dy 



7/2 - A2 1 + y . /f 



1 + eT 



2(1 -A2) 
1 



In 



(l + VT^y^T^ - Act) (l + VT^^/l^+ Act 



A2 



In 



1 + VT^y^T^ + Xer) {er - A) 



'T4 



1 

Av^r^ (^l + yr^yr^^ _ Ae^) (e^ + A) 
y^ dy 



(16) 



2-A2 



arctanh A 



l-e2. 



+ 



1-A2 ' 2(1 -A2) e2^-A2 



(17) 



To compute frp — tx we have to subtract huge numbers that are almost identical. Therefore 
we develop the relevant differences of the four integrals separately. Setting x := a' /a ~ 
Tg/ro we have: 



An 

GM 

/S.j'2 
At4 



-(l/a;-l) 



arctanh (A/ct) 



Tl 



i(l-x^)^-i(l/x^-l)^ "^^^""y/^^^ 
A 



^^B^)rT2 - VWo)It2 ~ In X, 



(l/x-l) 



arctanh(A/eT) 
A 



+ Inx - - 1)6 



arctanh(A/eT) 
A ^ 



/BK)/^4 - V5(ro)/T4 ~ In 



(18) 

(19) 

(20) 
(21) 
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The huge number — l)arctanh(A/eT)/A cancels when we add the four terms: 
t'r^-tx ~ An + GMAt2 + GMAts + 2 GMAn 
~ CM 



1(1 - x^) y - - 1)5 - 2 Inx 



A 



Finally we have the time delay to leading order in 6, e. and A, 



At ~ Jl- Ar^/3GM 



1 _ 

2 1 



1+^ 



- Ar2/3 GM 

GM ( , /a , /a 

— , — arctanhA/ — + arctanhA/ — ro 

V^r2 \^ V 3 V 3 ^ 



+41n — 
a' 



:22i 



(23) 



4 SDSS J1004+4112 

Consider the lensing cluster of SDSS J1004+4112 and the quasar as source [HIE]. As we 
have at least 4 images, the cluster cannot be spherically symmetric. We will close our eyes 
to the images with small a and consider only the images C and D with a = 10" ± 10 % 
and a' = 5" ± 10 %. We argue that they are less senitive to a non-spherical inner 
struture of the cluster. The mass of the cluster is M = (1 ± 0.2) ■ 10^^ kg. The cluster 
has a redshift of = 0.68. A numerical integration of the Cold Dark Matter model with 
A = 1.5 ■ 10~^2 m~2 yields an area distance = 7.0 ■ 10^^ m as seen from Earth from. 
This area distance coincides with the coordinate distance r^. For the quasar we have 
Zs = 1.734 yielding the area distance ds = 13.7 ■ 10^^ m. Due to the magnification of the 
quasar by the cluster, the translation from its area distance to its coordinate distance is 
ambiguous. We will use two Ansatze, 



d, 



rr + rs 
1 - Ar|/3 



or 



ds = rT + rs- 



(24) 



The first one gives rs = 5.3- 10^^ m, A = 2.2- 10"^^ m"^ from equation (fTOl) and At = 20.1 
years. The second gives rs = 6.7 • 10^^ m, A = 2.7 ■ 10"^^ m"^ and At = 17.1 years. The 
errors on these quantities coming from the errors in M, a and a' can be read from tables 
1 and 2, 



A = 
At = 



(2.0 ± 1.4) ■ 10 
(20.2 ± 7.5) years 



52 2 



or 
or 



(2.4 ± 1.5) ■ 10 
(18.1 ±6.1) years. 



52 m-2, 



(25) 
(26) 



It is encouraging to note that both constraints fl25|) on the cosmological constant are 
compatible with the present observational bounds, A = (1.5 ± 0.7) ■ 10"^^ m"^, whose 
central value we have used without variation to translate redshifts into area distances. 
Also both constraints (!26l) on the time delay are compatible with the recent observational 
lower bound [1] At > 5.7 years. 
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M ± 20% 


+ 










+ 


+ 


+ 




+ 


+ 






+ 


+ 






a' ± 10% 


+ 


+ 


+ 










+ 


rs [10^5 m] 


5.3 


6.1 


5.6 


5.3 


5.6 


5.0 


4.8 


5.0 


A[10-52 m-2] 


2.2 


0.7 


1.5 


2.2 


1.5 


2.8 


3.4 


2.8 


Ar [years] 


24.1 


20.3 


12.8 


16.1 


23.6 


27.7 


18.7 


15.1 



Table 1: Values of from the first of equations (1241) . A and Ar at the corners of the 
error box in M, a and a' . Minimal and maximal values are bold face. 



M ± 20% 


+ 










+ 


+ 


+ 


a ± 10% 


+ 


+ 






+ 


+ 






a' ± 10% 


+ 


+ 


+ 










+ 


A[10-52 m-2] 


2.7 


0.9 


1.9 


2.7 


1.9 


2.3 


3.9 


3.3 


Ar [years] 


21.8 


19.6 


12.0 


14.6 


22.0 


24.2 


16.0 


13.3 



Table 2: Values of A and Ar for rg = ds — tt = 6.7 ■ 10^^ m at the corners of the error 
box in M, a and a' . Minimal and maximal values are bold face. 

5 Conclusions 

Our computation rests on three shaky idealisations: 

• the sphericity of the cluster, 

• ignoring the velocity of observer and source with respect to the lens, 

• ignoring all masses except for the cluster. 

The first raises the good old question whether a spherical cow can be useful. To go beyond 
it implies to take into account multiple scattering of the photon off sub-constituents of 
the cluster. Calculating angles, time delays and number of images from the density profile 
of the cluster and its dark matter halo is a well developed art [lO]. Applying it to SDSS 
J1004+4112, Kawano & Oguri [11] predict a time delay between images C and D of up 
to 10 years. However including a positive cosmological constant in their analysis is not 
straight forward. 

Khriplovich & Pomeransky [12] point out that if the earth is taken comoving with 
respect to the exponentially expanding de Sitter space then the effect of the cosmological 
constant on the deflection cancels. 

Going beyond the third idealisation necessitates an interpolating solution, which em- 
beds many static, curved Kottler solutions into the ambient expanding, fiat Friedmann 
solution. This is a long standing problem, to which already Einstein & Straus [131 E] 
have contributed an unstable [15] solution. A first qualitative assessment based on this 
solution is given in section III of reference [1^: the other clusters in the universe weaken 
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the effect of tlie cosmological constant and increase upper bounds on the cosmological 
constant from certain lenses by two orders of magnitude. 

It is safe to conclude: we all agree that lensing and time delay at cosmological scale 
is beautiful physics and that our theoretical understanding of it is incomplete. For the 
older physicists among us however, a time delay of 10 or 20 years is bad news. 

Note added: Mustapha Ishak [17j has extended the analysis of how the embedding of 
Kottler's solution into Friedmann's solution modifies light deflection to include the time 
delay. We thank him for having kindly sent us his work. 
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